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1  Introduction 

Several  major  problems,  namely,  uncertain  surface  forcing  function,  unknown 
open  boundary  conditions  (OBC),  and  pressure  gradient  error  using  the  er- 
coordinate,  affect  the  accuracy  of  coastal  ocean  prediction.  At  open  lateral 
boundaries  where  the  numerical  grid  ends,  the  fluid  motion  should  be  unre¬ 
stricted.  Ideal  open  boundaries  are  transparent  to  motions.  The  most  popu¬ 
lar  and  successful  scheme  is  the  adjoint  method.  The  disadvantages  that  may 
restrict  its  use  are  ocean-model  dependency  and  difficulty  in  deriving  the  ad¬ 
joint  equation  when  the  model  contains  rapid  (discontinuous)  processes,  such 
as  change  of  ocean  mixed  layer  from  entrainment  to  shallowing  regime.  De¬ 
velopment  of  a  ocean-model  independent  algorithm  for  determining  the  OBC 
becomes  urgent. 

Reduction  of  horizontal  pressure  gradient  error  is  another  key  issue  of  using 
^-coordinate  ocean  models,  especially  of  using  coastal  models.  The  error  is 
caused  by  the  splitting  of  the  horizontal  pressure  gradient  term  into  two  parts 
and  the  subsequent  incomplete  cancellation  of  the  truncation  errors  of  those 
parts.  As  advances  of  the  computer  technology,  use  of  highly  accurate  schemes 
for  ocean  models  becomes  feasible. 


2  Uncertainty  of  Surface  Forcing 

To  investigate  the  uncertainty  of  surface  wind  forcing  and  its  effect  on  the  coastal 
prediction,  the  Princeton  Ocean  Model  (POM,  Blumberg  and  Mellor,  1987) 
was  used  with  20  km  horizontal  resolution  and  23  sigma  levels  conforming  to  a 
realistic  bottom  topography  during  the  life  time  of  tropical  cyclone  Ernie  1996 
over  the  South  China  Sea  (SCS).  A  study  (Chu,  et  al,  1999)  shows  that  the 
root-mean-square  (RMS)  difference  of  each  component  (zonal  or  latitudinal) 
between  the  two  wind  data  (NCEP  and  NSCAT)  over  the  whole  SCS  during 
November  1996  fluctuated  between  2.7  m/s  to  6.5  m/s.  The  uncertainty  of  the 
whole  SCS  response  to  the  two  wind  data  sets  were  4.4  cm  for  surface  elevation, 
0.16  m/s  for  surface  current  velocity,  and  0.5°C  for  near-surface  temperature, 
respectively. 
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3  Jacobian  Matrix  Method  for  Determining  OBCs 

Improvement  of  coastal  prediction  largely  depends  on  determination  of  lateral 
OBCs,  vector,  B  =  (61,62,  ...,6n).  The  observation  forms  an  m-dimensional 
vector  (observation  vector)  O  —(0i,02,  —Om),  located  at  the  interior.  If  B 
is  given,  we  can  solve  the  dynamic  system  and  obtain  the  solution  S.  At  the 
observational  points,  the  solutions  form  a  solution  vector  S  =  (Si, S2, 
which  depends  on  B  (Fig.  1). 

It  is  reasonable  to  determine  B  with  the  given  O  by  minimize  the  RMS  error 


1  m 

J=  -^(Sj-Oj)2. 

\mU 


(1) 


which  leads  to  a  set  of  n  equations  implicitly  solvable  for  61,62,  ...&„, 


~  Oj)Rij  =  0,  *  =  1,2,...,  1 

j=  1 


(2) 


where 


Rij  =  *  =  1)2, ...,  n;  j  =  1, 2,  ...m.  (3) 

are  components  ofanxm  Jacobian  matrix  R  =  {Rij}- 

From  a  first  guess  boundary  vector  B*,  a  solution  vector  S*  is  obtained  by 
solving  the  numerical  ocean  model.  The  RMS  between  S*  and  O  might  not  be 
minimal.  We  update  the  boundary  parameter  vector  components  by  increments 
{<16,  |  i  =  1, 2 . n} ,  and  therefore  components  of  the  solution  vector  become 


Sj  =  S*  +  y^  RijSbi  4-  high  order  terms  (4) 

t=l 

Substituting  (4)  into  (2)  and  neglecting  higher  order  terms  leads  to  a  set  of  n 
linear  algebraic  equations  for  {dbi}  , 

n 

Y^PuSbi=di,  *  =  1,2,  ...,n  (5) 

i= 1 

where 

m  m 

Pji  =  Rij  Rij ,  dj  —  ^  Rjj  {Oj  —  Sj );  *  =  1,2, ...,  n;  1  =  1,2,  ...n. 

j— 1  j=i 

(6) 

Both  Oj  and  Sj  are  known  quantities.  Therefore,  the  linear  algebraic  equations 
(5)  have  definite  solutions  when  the  Jacobian  matrix  {Rij}  is  determined  and 


det  {Pu}  i1  0 


(7) 
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This  method  was  verified  by  a  flat  bay  centered  at  35°  N  and  bounded  by 
three  rigid  boundaries  (Fig.  1).  This  bay  expands  1000  km  in  both  the  north- 
south  and  east- west  directions.  The  northern,  southern,  and  western  boundaries 
are  rigid,  and  the  eastern  boundary  is  open.  Using  the  optimization  method, 
the  temporally  varying  OBC,  B (t),  is  determined.  After  10  day’s  of  integration, 
the  magnitude  of  relative  error 


p(0)  -  Z  I  Oj  ~  Sjl 

■  mi 


(8) 


is  on  the  order  of  10-4-10-5  (Fig.  2),  which  is  almost  in  the  noise  level.  The 
Jacobian  matrix  method  performs  well  even  when  random  noises  axe  added  to 
the  ’observational’  points.  This  indicates  that  we  can  use  real-time  data  to 
invert  for  the  unknown  open  boundary  values. 


4  High-Order  Difference  Schemes 

4.1  A  Hidden  Problem 

Improvement  of  the  prediction  also  partially  depends  on  the  selection  of  the  dis¬ 
cretization  schemes.  Most  coastal  models  use  second-order  difference  schemes 
(such  as  second-order  staggered  C-grid  scheme)  to  approximate  first-order  deriva¬ 
tive 


( <?P\  „  Pi+ 1/2  Pi— 1/2  _}_(  d3p\  2 
\dx )i~  A  24\dx3)i  ’ 


(9) 


where  p,  A  represent  pressure  and  grid  spacing.  This  scheme  uses  the  local 
Lagrangian  Polynomials  whose  derivatives  are  discontinuous. 


4.2  Combined  Compact  Scheme 

Recently,  Chu  and  Fan  (1997,  1998,  1999,  2000)  proposed  a  new  three-point 
combined  compact  difference  (CCD)  scheme, 


—  (/*+i  -  2 fi  +  fi- 1) 


(10) 
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to  compute  f[  ,  f by  means  of  the  values  and  derivatives  at  the  two 
neighboring  points.  Moving  from  the  one  boundary  to  the  other,  CCD  forms  a 
global  algorithm  to  compute  various  derivatives  at  all  grid  points,  and  guaran¬ 
tees  continuity  of  all  derivatives  at  each  grid  point. 

4.3  Seamount  Test  Case 

4.3.1  Model  Description 

Suppose  a  seamount  (Fig.  3)  located  inside  a  periodic  /-plane  (  /o  =  10~4s-1) 
channel  with  two  solid,  free-slip  boundaries  along  constant  y.  Unforced  flow 
over  seamount  in  the  presence  of  resting,  level  isopycnals  is  an  idea  test  case 
for  the  assessment  of  pressure  gradient  errors  in  simulating  stratified  flow  over 
topography.  The  flow  is  assumed  to  be  reentrant  (periodic)  in  the  along  channel 
coordinate  (i.e.,  x-axis).  We  use  this  seamount  case  of  the  Semi-spectral  Prim¬ 
itive  Equation  Model  (SPEM)  version  3.9  (Haidvogel  et  al.,  1991)  to  test  the 
new  difference  scheme. 

4.4  Temporal  Variations  of  Peak  Error  Velocity 

Owing  to  a  very  large  number  of  calculations  performed,  we  discuss  the  results 
exclusively  in  terms  of  the  maximum  absolute  value  the  spurious  velocity  (called 
peak  error  velocity)  generated  by  the  pressure  gradient  errors.  Fig.  4  shows  the 
time  evolution  of  the  peak  error  velocity  for  the  first  20  days  of  integration 
with  the  second-,  fourth-,  and  sixth-order  ordinary  schemes.  The  peak  error 
velocity  fluctuates  rapidly  during  the  first  few  days  integration.  After  the  5 
days  of  integration,  the  peak  error  velocity  show  the  decaying  inertial  oscillation 
superimposed  into  asymptotic  values.  The  asymptotic  value  is  around  0.19  cm/s 
for  the  ordinary  scheme  and  0.15  cm/s  for  the  compact  scheme.  For  the  sixth 
order  difference  the  asymptotic  value  is  near  0.04  cm/s  for  the  ordinary  scheme 
and  0.02  cm/s  for  the  compact  scheme. 

5  Conclusions 

(1)  The  surface  forcing  function  contains  uncertainty.  The  difference  between 
commonly  used  NSC  AT  and  NCEP  surface  wind  data  is  not  negligible.  The 
response  of  the  South  China  Sea  to  the  uncertain  surface  forcing  is  also  evident. 
Therefore,  it  is  quite  urgent  to  study  the  model  sensitivity  to  surface  boundary 
conditions. 

(2)  The  Jacobian  matrix  method  provides  a  useful  scheme  to  obtain  unknown 
open  boundary  values  from  known  interior  values.  The  optimization  method 
performs  well  even  when  random  noises  are  added  to  the  ’observational’  points. 
This  indicates  that  we  can  use  real-time  data  to  invert  for  the  unknown  open 
boundary  values. 

(3)  The  er-coordinate,  pressure  gradient  error  depends  on  the  choice  of  differ¬ 
ence  schemes.  Fourier  analysis  shows  that  the  fourth-order  scheme  may  reduce 
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the  truncation  errors  by  1-2  order  of  magnitude  compared  to  the  second-order 
scheme,  and  the  sixth-order  scheme  may  reduce  the  truncation  errors  further 
by  1-2  order  of  magnitude  compared  to  the  fourth-order  scheme.  Within  the 
same  order  of  the  difference  the  combined  compact  scheme  leads  to  a  minimum 
truncation  error.  The  compact  scheme  may  reduce  near  55%  error,  and  the 
combined  compact  scheme  may  reduce  near  84%  error  compared  to  ordinary 
sixth-order  difference  scheme. 
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Fig.  1 .  Determination  of  open  boundary  condition 
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Fig.  3.  Seamount  geometry  Fig.  4.  Peak  error  velocity  for  different  order  schemes 


